Spin-Orbit Coupling and Symmetry of the Order Parameter in Strontium Ruthenate 
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Determination of the orbital symmetry of a state in spin triplet Sr2Ru04 superconductor is a 
challenge of considerable importance. Most of the experiments show that the chiral state of the 
z(k x ± iky) type is realized and remains stable on lowering the temperature. Here we have studied 
the stability of various superconducting states of S^RuCU in the presence of spin-orbit coupling. 
Numerically we found that the chiral state is never the minimum energy. Alone among the five states 
studied it has (L.S) = and is therefore not affected to linear order in the coupling parameter A. We 
found that stability of the chiral state requires spin dependent pairing interactions. This imposes 
strong constraint on the pairing mechanism. 



Despite 10 years of intensive research [l|, the mech- 
anism of superconductivity in strontium ruthenate 
Sr2Ru04 is still not clarified [HQ- A necessary first 
step in identifying the effective interactions responsible 
for the superconductivity is to fully identify the sym- 
metry of the pairing state and to explain the various 
relevant experimental properties below T c . The cur- 
rent consensus is that most experiments point towards 
a p-wave spin triplet, nodeless chiral state of the type 
d(k) = z(sin& x + isinfcj,), but with additional lines of 
gap nodes somewhere on the Fermi surface 0, Q . In the 
Zhitomirsky and Rice model 5] the lines of gap nodes 
arise because of an inter-band proximity effect which 
couples of the Cooper pairs between the "active" and 
"passive" sheets of the Fermi surface. In our alternative 
proposal 0, we introduced an effective interlayer cou- 
pling model, which also leads to a chiral symmetry state 
with lines of gap nodes. This model has a fully gaped 
chiral state on the 7 Fermi surface sheet and horizon- 
tal lines of nodes at k z = ±ir/c on the /3 sheet of the 
Fermi surface. Moreover, we have shown that it is in 
good agreement with the experimental measurements of 
specific heat, penetration depth and thermal conductiv- 
ity below T c . 

In this letter we address the crucial question of the 
stability of the chiral symmetry state in the presence of 
spin-orbit coupling. For p-wave spin-triplet superfluids 
a large number of different pairing states are possible, 
and in general many different states are degenerate in 
the weak coupling limit. In the well known case of su- 
perfluid 3 He, weak coupling theory favors the fully gaped 
Balian-Werthammer (B) phase relative to the Anderson- 
Brinkman-Morrel (A) phase which has point gap nodes 
||. In this case the A phase is stabilized (in zero mag- 
netic field) only by a spin-fluctuation feedback effect in 
which the strength of the effective pairing interaction it- 
self is modified by the pairing symmetry. The broken 
symmetry ground state effectively leads to a magnetic 
anisotropy in the spin fluctuation spectrum, which in 
turn favors the chiral solution. The analogous question 



in Sr2Ru04 is to ask whether the experimentally deter- 
mined chiral state is stable simply because of magnetic 
anisotropy arising from spin-orbit interactions in the nor- 
mal state, or whether it is also necessary to invoke a di- 
rectly spin-dependent pairing interaction. Below we ar- 
gue that this latter case applies, and therefore the very 
existence of the chiral pairing state provides evidence for 
a magnetic contribution to the pai ring mechanism. 

We follow earlier calculations |3, [lOj, |llj, |l2| where the 
role of spin-orbit coupling in the superconducting and 
normal states, respectively have been investigated. Our 
results are qualitatively consistent with the supercon- 
ducting state calculations of Ng and Sigrist 0, [ul Hj| , 
but, importantly, differ in conclusion. We show below 
that single-particle spin-orbit coupling alone is unlikely 
to stabilize the chiral state. The chiral state is stabi- 
lized only after allowing for a spin-dependent effective 
interaction. Such an effective interaction is consistent 
with the implication of spin-orbit interaction for trans- 
verse and longitudinal spin fluctuations, in the work of 
Eremin, Manske and Bennemann |l^ and Kuwabara and 
Ogata 13]. Therefore we suggest that the experimental 
observation of a chiral pairing state is conclusive proof 
of a spin-dependent pairing interaction. Whilst this con- 
clusion strongly supports the spin fluctuation mechanism 
of pairing, it could also be consistent with other spin- 
dependent interaction models, such as Hund's rule cou- 
pling Q. 

The effective pairing Hamiltonian we consider is a sim- 
ple three-band, attractive U, Hubbard model: 



H = ((e m - n)SijS mm ' -tmm'(ij))c. 

ijmm' ,cr 



"imijCjni'a 



' r, 5^ ^mrrJ (^j)^ima^jm' p^jm' '■yCimS , (1) 



where m and ml refer to the three Ruthenium t2 g orbitals 
a = xz, b = yz and c — xy and i and j label the sites of a 
body centered tetragonal lattice. The hopping integrals 
tmm' (ij) and site energies e m were fitted to reproduce the 
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FIG. 1: Critical temperature for solutions of various symme- 
tries denoted by symbols (a)-(e) defined in Eq. (5). 



experimentally determined Fermi surface |l5l Il6| . In our 
previous papers 0, we chose the simplest set of attrac- 
tive pairing interactions U mm /(ij) which lead to results 
consistent with experiment. However, more generally, for 
a spin-dependent pairing interaction the effective Hub- 
bard parameters Um^J S (ij) are spin as well as orbital 
dependent. 

In the presence of spin-orbit coupling the original 
Hamiltonian (Eq. f ) should be supplemented by the term 
AL.S for each Ru atomic site. This can be expressed 
in the form 



H 



so 



. A 



mm _k _+ 



(2) 



i.erer' mm' 



where a£ a , , k = x,y, z, are the Pauli matrices, e Kmm de- 
notes the completely antisymmetric tensor, and the sign 
convention implies that here the Ru orbital indices must 
be ordered asm = {yz, zx, xy) or (&, a, c) in our notation 


In the Hartree-Fock-Gorkov approximation the gap 
equation takes the simplified form 

Aw (y ) = E>w (<i )Xw (v) • (3) 

where the pairing amplitude, or order parameter, 

Xmm'fe)> i s defined by the usual relations [ijfl 

Xl-(ii) = ^< Bff < v (l-2/(^)) for a = a' 
1 



Xmm' 



\ zm<J w jm'<r' 



(4) 



-vZ„'U , < m , <T )(l-2f(E»)) for <7^a'. 



Note that the full spin and orbital dependent effec- 
tive interaction U^J S (ij) from Eq. is reduced to 
the simpler set of parameters U^ m ,(ij) for which there 
are only two distinct cases: llVi = and 
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FIG. 2: Condensation free energy for solutions of various 
symmetries defined in Eq. 0. For A > the most stable 
state is (b) d(k) = (X, Y, 0) , while for A < the state (e) 
d(k) = (Y,X, 0)). The chiral state (a) is never stable for 
nonzero spin-orbit coupling. The inset shows the condensa- 
tion energy for small A calculated for the states (b)-(e) with 
respect to the state (a). 



Uli m '(ij) = U^ m ,(ij). Other more general spin depen- 
dent terms, such as spin-flipping terms, might occur, but 
these would not have any effect on the states we consider. 

The group theory of pairing symmetry of triplet (odd- 
parity) superconductivity in tetragonal crystals has been 
discussed by Volovik and Gorkov |l8( , and is reviewed in 
Annett 0] an d Mineev and Samhokin |2o|. In this case 
there are five distinct symmetry gap functions which are 
relevant. Transforming to k-space and using the stan- 
dard representation A£ CT = ia y a.d. k , where A and d are 
matrices in orbital indices, they have the characteristic 
order parameters 
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corresponding to the E u , A\ u , A<2 U , B\ u and i?2« pairing 
symmetry states [l^. Here X(k) and T(k) are pairs of 
basis functions transforming as k x and k y , respectively, 
under crystal point group symmetries. For in-plane in- 
teractions X = sin k x and for inter-plane interactions 
X = sin (k x /2) cos (k y /2) cos (k z c/2) in the body cen- 
tered tetragonal crystal Sr 2 Ru04. For brevity we shall 
refer to these five pairing states as (a)-(e) below. 

If the pairing interaction were spin- independent, then 
the interaction parameters would obey 



U^(ij) = U^ m ,(ij). (6) 
With this assumption we solved the full self-consistent 
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gap equations (Eqs. in presence of a weak spin- 

orbit coupling parameter A. The band structure and the 
interaction par ameters were chosen as in our previous 
works |g, |21| . The results for the critical temperature 
T c are shown in Fig. n an d the condensation free energies 
(F s — F n ) at T = are shown in Fig. |21 Both figures 
show clearly that all five pairing states become exactly 
degenerate at A = 0, as expected. For positive values 
of A the states split in such a way that the (b) and (c) 
states have the highest T c and the greatest (negative) 
condensation energy. Conversely, for negative values of 
A states (d) and (e) have the highest T c and the greatest 
condensation energy. Surprisingly, the chiral state, (a), 
is never stable compared to the others. 

This counter- intuitive result deserves some comment. 
Firstly, one should note that experiments point unequiv- 
ocally to a chiral solution, since only this one would have 
the constant in-plane spin susceptibilit y as measured in 
both NMR 22] and neutron scattering |23j experiments. 
Secondly, Ng and Sigrist 0, 0, have argued on the 
basis of a very similar model to ours that spin-orbit in- 
teraction does stabilizes the chiral phase. 

Although we draw different conclusions, on closer in- 
spection, our results are indeed very similar to theirs. 
They used a three-band model characterized by two main 
pairing parameters V\ = Vi and V3 . They found that the 
chiral state is favored for positive A only in the region 
V3/V2 < 1. For V3/V2 > 1 they found that states (b) and 
(c) (III and II in their notation) were favored, and that 
these two states were nearly degenerate being separated 
only by an additional inter-band proximity effect cou- 
pling term. Our results are qualitatively consistent with 
this finding, if we interpret our effective Hubbard inter- 
action terms 0, U± and U\\ as giving a corresponding 
value of V3/V2 « 2U±/U\\ = 2.4. In our model we have 
no inter-band proximity effect coupling, and hence states 
(b) and (c) are essentially degenerate. 

The physical values of the spin-orbit coupling param- 
eter are positive, and believed to be of order A = 40 — 
80meV, corresponding to X/t ~ 0.5 — 1 in our units. We 
have also carried out fully self-consistent calculations at 
A = 0.5, and found the same ordering of the ground states 
as shown in Fig. [21 But for such large values of A we ob- 
served that the Fermi surface shape begins to change, 
and so it was necessary to refit the Hubbard model hop- 
ping parameters to the experimental Fermi surface 
However this refitting of the Fermi surface did not change 
the results qualitatively from those expected from Fig. |3 
and so we conclude that the precise Fermi surface shape 
is not essential in determining the relative stability of the 
different pairing states. 

In order to understand more fully the physical origin 
of the spin-orbit splitting shown in Figs. we consider 
small perturbations about the point A = 0. Here the 



change in condensation energy to the linear order is 
«5F = ^(L.S) = ^(m|4|m>TT ro , 

+ \ J2 H^* - i£,/l"*>w (7) 

m,m' 
m,m' 

where n^ m , = (cf ma Ci m i „/) are on-site density matrices 
for the Ru d orbitals, in the absence of spin-orbit cou- 
pling. In fact, in our model, the 7 band associated with 
the d xy orbital is strongly decoupled from the d xz and 
d yz orbitals, and so the only significant off-diagonal den- 
sity matrix elements are and n%£ . Thus we find the 
spin-orbit correction 

dF = ~i(n%-nit-nli+nti). (8) 

For the chiral (a) state n^J b = nil and so the spin-orbit 
coupling is zero to linear order in A. On the other hand 
for (b)-(e) nil = ~~ n al ana - so ^ ne spin-orbit energy is of 
order A. This explains the small A behavior seen in our 
numerical calculations, and shown in the inset to Fig. 
It is interesting to note that for all the states (a)-(e) the 
expectation value of the spin (S) = for each orbital 
to. On the other hand the z component of the orbital 
angular momentum is 

+ n l a l b -(41 +nH), (9) 

which is non-zero only for the chiral state. Therefore 
(L.S) is non-zero for four states (b)-(e) which have (S) = 
(L) = while the (a) state has (L.S) = 0, (L z ) ^ 0. 

Given this surprising result, how we can understand 
the stability of the chiral state, which is required for un- 
derstanding the susceptibility experiments? Clearly, if 
the effective pairing interaction is caused, even in part, by 
spin fluctuations, then in general we would expect that 
^mm' (*■?') ^ ^i B <(ij)i corresponding to the two Fermi 
liquid parameters F^ , F^ |8j. In the spin-fluctuation 
theory of Kuwabara and Ogata 0] one can see that 
the parallel spin and opposite spin effective interactions 
are derived from transverse or longitudinal spin fluctua- 
tions, respectively, and hence they become different when 
spin-orbit interactions are present. The role of spin- 
orbit coupling in the transverse, x ± (q), and longitudinal, 
X 2 (q), spin susceptibilities has been exp lored more fully 
by Eremin, Manske and Benneman n 1121 and their theory 
agrees well with the results of NMR |22( and neutron scat- 
tering experiments |24l |. Therefore spin-orbit coupling 
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FIG. 3: The same as in Fig. |21 but here we assumed that 
interactions for 'a' solution U are slightly (by 1%) larger than 
for other solutions U: U' /U = 1.01. 



mal state of Sr2Ru04. They found that spin-orbit in- 
teraction led to a magnetic anisotropy between \ z (q> w ) 
and x ± (q, w). The well known argument for 3 He would 
then suggest that these susceptibilities lead to different 
effective pairing interactions between spin-triplet, S = 1, 
quasiparticle pairs in the m = or m = ±1 chan- 
nels. As it turns out a small increase of the interaction, 
U' /U = 1.01, stabilizes the (a) state. 
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Linkage Grant No. 979446 and the INTAS grant No. 
01-654. GL would like to thank Max Planck Institute for 
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ity. The authors are grateful D. Manske and I. Eremin 
for discussions. 
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FIG. 4: Phase diagram for different interaction ratio U' /U 
and a spin-orbit coupling A. 



has two quite distinct effects on a spin triplet supercon- 
ductor. Firstly, the spin-orbit term, Eq. J5J, directly 
enters the gap equation leading to small changes in the 
single-particle band structure and density of states, and 
secondly, it leads to spin anisotropy in the effective pair- 
ing interaction. Both of these effects must be included 
in order to establish the relative stability of the different 
possible pairing symmetry states |2fij . 

We have found that the chiral state becomes stable 
when we increased the interaction = U' by 1% with 
respect to W* = JJ. These results in Fig. [31 The phase 
diagram in the U'/U vs. X/t space is presented in Fig. 0] 
It shows that to stabilize the chiral state in the presence 
of sizeable spin-orbit coupling requires increase of U' by 
less than 1%. 

It follows from these results, that the pairing inter- 
action itself must directly favor the chiral symmetry 
state and would be an analogue of the well known spin- 
fluctuation feedback mechanism which stabilizes the A- 
phase in superfluid 3 He 8]. Eremin, Manske and Bcnnc- 
mann calculated the spin susceptibility in the nor- 
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